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EMPIRICAL MODELS FOR USE IN DESIGNING 
DECOMPRESSION PROCEDURES FOR SPACE OPERATIONS 


ABSTRACT 

Empirical models for predicting the incidence of Type I altitude decompression sickness (DCS) and 
venous gas emboli (VGE) during space extravehicular activity (EVA), and for use in designing safe 
denitrogenation-decompression procedures are developed. The parameters of the models are esti- 
mated using DCS and VGE incidence data from 26 NASA and USAF manned altitude chamber 
decompression tests involving 607 male and female subject tests. 

These models, and procedures for their use, consist of (1) an exponential relaxation model and pro- 
cedure for computing tissue nitrogen partial pressure resulting from a specified prebreathing and 
stepped decompression sequence; (2) a formula for calculating Tissue Ratio (TR), a tissue decom- 
pression stress index; (3) linear and Hill equation models for predicting the total incidence of VGE 
and DCS attendant with a particular TR; (4) graphs of cumulative DCS and VGE incidence (risk) 
versus EVA exposure time at any specified TR; and (5) two equations for calculating the average 
delay period (latency time) for the initial detection of VGE or indication of Type I DCS in a group 
after a specific denitrogenation-decompression procedure. Several examples of realistic EVA 
preparations are provided to demonstrate the use of the predictive models and interpretation of 
the results. 


INTRODUCTION 

This report documents empirical models based on NASA and USAF altitude chamber experiments 
for predicting the incidence (risk) of Type I altitude decompression sickness (DCS) and venous gas 
emboli (VGE) during space EVA operations, and for use in designing safe denitrogenation-decom- 
pression procedures consistent with an acceptable level of DCS and VGE risk. The models enable 
predictions, with a quantitative level of confidence and accuracy of the incidence of DCS and VGE 
to be expected within a group as a result of a particular prebreathe-decompression procedure. 

Earlier guidelines for predicting DCS were based on data collected from a literature survey on past 
altitude DCS research. No data on VGE are available from those reports since the instrumentation 
for in vivo gas phase detection had not been developed. In addition, the type, intensity, duration 
of exercise, and altitude exposure time in those studies were not consistent with NASA's current 
EVA procedures. 

The mathematical models described consist of (1) a time-dependent exponential model for computing 
tissue nitrogen partial pressure under stepped prebreathe and decompression conditions, (2) a form- 
ula for computing tissue nitrogen decompression stress or Tissue Ratio (TR), and (3) empirical linear 
and Hill equation models for predicting the incidence of Type I DCS, VGE, and Type I DCS Grade 3 
(symptoms that forced early test termination or that reappeared after the test). 
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The predictive models are obtained by mathematically fitting the linear and Hill equations to the 
paired TR and % DCS or % VGE data obtained from 26 unique NASA and USAF manned altitude 
chamber tests that simulated EVA work profiles. In total, 607 male and female exposures were 
performed from 1982 to 1986 under controlled conditions at the Johnson Space Center, Houston, 
Texas, and Brooks Air Force Base, San Antonio, Texas. 


DEVELOPMENT OF TISSUE DECOMPRESSION STRESS (TISSUE RATIO) MODEL 


Nitrogen Uptake and Washout Exponential Model 

All body tissues in equilibrium with ambient nitrogen will respond to increased or decreased nitrogen 
partial pressure in the breathing medium by absorbing or eliminating the gas that is dissolved in 
tissues. The rate of elimination/absorption is determined by the rate constant of the tissue in ques- 
tion and the nitrogen partial pressure difference between the tissue and ambient nitrogen, i.e., 

p a _ P j (1) 


dP/dt = k 


where 


dP/dt = time rate of change of tissue nitrogen partial pressure, 

P = tissue nitrogen partial pressure at any time t, 

P a = ambient nitrogen partial pressure, 

k = tissue nitrogen partial pressure rate constant. 

The rate constant (k) of a given tissue is physiologically determined by the blood perfusion rate 
through the tissue, the diffusion rate of the inert gas through the tissue, and the solubility char- 
acteristics of the gas in the specific tissue. 

The tissue rate constant (k) is related to the tissue nitrogen half-time (t 1/2 ) through the equation. 


k = (In 2)/ 1 = .693/t 1/9 

1/2 1/2 


where 


ti/g = tissue nitrogen partial pressure half-time 
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The body can be characterized by a spectrum of theoretical tissue types that represent "fast” and 
"slow” tissues. A theoretical tissue type is characterized by the tissue half-time (t 1/2 ): the time re- 
quired for a tissue to respond to a change in ambient nitrogen partial pressure by giving off or 
absorbing nitrogen until the initial difference between the tissue nitrogen partial pressure and 
ambient nitrogen partial pressure is reduced by one-half. 

Since t 1/2 characterizes the nitrogen perfusion-diffusion properties of the tissue, k is also a charac- 
teristic of the tissue-nitrogen relationship. It can be seen from equation 1 that, for a given (P - P) 
difference, a tissue having a greater k- value will absorb or eliminate nitrogen more rapidly than 
one with a lesser k-value. A "fast” tissue such as blood may have a half-time of only a few minutes 
while a "slow" tissue half-time may be several hours. 

Following a rapid change in nitrogen partial pressure in the breathing medium, the nitrogen par- 
tial pressure which is reached in a designated "tissue type" after a specific time is described by equa- 
tion 2, the solution of equation 1. This equation describes either nitrogen uptake or elimination in 
a specific theoretical tissue type depending on the initial tissue nitrogen partial pressure and the 
partial pressure of nitrogen in the ambient breathing mixture. 

p .= p . + ( p .- p .)( i -' -kt ) <a 


where 

P t = the nitrogen partial pressure in the tissue after exposure for t minutes, 

P Q = initial tissue nitrogen partial pressure, 

P a = ambient nitrogen partial pressure in breathing medium, 

e = base of natural logarithm, 

t = exposure time in minutes, and 

k = tissue nitrogen partial pressure rate constant. 

The exponential decay (relaxation) model (eq.2) is generally used to describe the final inert gas 
(in this case nitrogen) partial pressure within a specific theoretical body tissue type after a par- 
ticular prebreathe procedure. Examples of its use are found in references 1 and 2. 


Nitrogen Super saturation and the Tissue Ratio Concept 

Tissue nitrogen supersaturation can be defined in two ways. The second definition is used in this 
report in constructing models that predict the incidence of DCS and VGE. 
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A person breathes 100% oxygen from a mask. This procedure is intended to reduce the total dissolved 
nitrogen in the body by eliminating the nitrogen concentration in the breathing mixture. This allows 
the body to "wash” nitrogen from its tissues through perfusion and gaseous diffusion processes. The 
tissues are saturated with nitrogen in relation to the breathing medium and nitrogen diffuses from 
the tissues to the lungs due to an imbalance in the nitrogen partial pressure. No DCS or VGE result 
from this tissue nitrogen saturation since ambient pressure is not reduced. The nitrogen remains 
dissolved in tissues and body fluids. 

Supersaturation exists when a system is in a metastable state and there is a potential for spontaneous 
change in phase of some constituent. This can occur when a stable equilibrium state is disturbed by a 
change in pressure, volume, temperature, or chemical constituent. Decompression (decrease in ambi- 
ent pressure) as experienced by divers or aviators produces a transient metastable state in the body 
tissues such that the rate of normal, random micronucleation events increase as well as the size of the 
gas nuclei. Two definitions of tissue nitrogen supersaturation are presented. 

(1) A tissue is supersaturated with dissolved nitrogen when the nitrogen partial pressure in the 
tissue is greater than ambient pressure. 

With this definition of tissue nitrogen supersaturation, ambient pressure is reduced to a point 
where it is less than the partial pressure of nitrogen in a particular tissue. The tissue is not said 
to be supersaturated until the ambient pressure is less than the nitrogen partial pressure in the 
tissue. This implies and assumes that no DCS symptoms or VGE production will result as long as 
the decrease in ambient pressure is not greater than the nitrogen partial pressure in a particular 
tissue, and that nitrogen is the only gas that needs to be considered in nucleation events. 

(2) A tissue is supersturated with dissolved nitrogen and has the potential to form a gas phase 
whenever ambient pressure is reduced. 

This implies that the potential to form gas emboli exist whenever ambient pressure is reduced 
regardless of the magnitude of pressure reduction. This definition does not assume that a finite 
degree of decompression can be tolerated. It assumes that all dissolved gases are potentially 
involved in nucleation events. 

It is intuitively evident that the greater the reduction in ambient pressure, the greater the 
potential for decompression sickness. It is possible to establish an index of decompression stress 
that relates tissue nitrogen partial pressure and ambient pressure changes by developing the 
concept of Tissue Ratio. 

The Tissue Ratio (TR) concept, a measure of tissue nitrogen supersaturation or decompression stress, 
was developed by Haldane 2,12 for divers and later modified for aviators. He observed that the 
human body could be rapidly decompressed from 2 to 1 atmosphere or from sea level (1 atmosphere) 
to approximately 20,000 feet (1/2 atmosphere) without symptoms of decompression sickness. Thus, 
it appeared to be safe to decompress aviators on a schedule which never allowed the total gas partial 
pressure in the tissues to exceed twice the ambient atmospheric pressure, a TR of 2.0. These observa- 
tions implied that the body could tolerate a degree of nitrogen supersaturation without overt mani- 
festations of decompression sickness. 
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Haldane's concept has been modified over the years as a more complete understanding of tissue gas 
exchange dynamics has been achieved and ultrasonic detection equipment has been developed to 
identify venous gas emboli (VGE) in situ . Using this technology, VGE have been detected in experi- 
mental human subjects decompressed to a pressure equivalent of 13,120 feet (Haldane's TR = 1.65) 
and in animals decompressed to only 6,560 feet (Haldane's TR = 1.26) 13 . Type I DCS symptoms 
have been reported from as little as 12,000 feet (Haldane’s TR = 1.57) 19 Thus, it is clear that 
Haldane TR's of less than 2.0 can evolve nitrogen bubbles (VGE), and produce signs and symptoms 
of decompression sickness. 

The calculation of the TR in equation 3 for a tissue of specific half-time uses the final tissue nitrogen 
partial pressure from one or more sequential applications of equation 2. 

Final Tissue pN 0 (abs) 

TR = (3) 

Ambient Pressure (abs) 

The TR calculation used in equation 3 is different from Haldane's TR in that the numerator is the par- 
tial pressure of dissolved nitrogen in the tissues prior to decompression rather than the total pressure. 
Initial equilibrium tissue pN 2 is taken to be ambient nitrogen pN 2 at 14.7 psia (P a = 11.6 psia). 

The resulting TR value for any theoretical tissue half-time type can be expressed as a single number 
and is unique for a particular prebreathe procedure regardless of the complexity of that procedure. 

This TR is associated (correlated) with the incidence of VGE and/or DCS occurring during the subse- 
quent decompression to establish a statistical relationship between a particular prebreathe procedure 
and the resulting VGE and DCS at that altitude. 

The potential to form a gas phase exists whenever ambient pressure is reduced. The second definition 
of tissue nitrogen supersaturation is therefore used in this report. Under equilibrium conditions at 
sea level pressure and gas composition the TR that represents no potential for DCS symptoms or VGE 
production is 11.6/14.7 = 0.78. Any reduction in ambient pressure will increase this TR if nitrogen 
washout procedures have not been initiated prior to the decompression. 


Example of Equation 2 and Tissue Ratio Application 

A simple example will provide a better understanding of how the exponential decay model and TR 
are used in the eventual development of the equations that predict DCS and VGE incidence in a pop- 
ulation that performs a specific prebreathe procedure. 

In this example, a person performs a 4.0 hour pure oxygen prebreathe while at 14.7 psia. Since no 
other prebreathe procedure preceded this event, all body tissues are initially in equilibrium with the 
ambient nitrogen partial pressure, P fl = 11.6 psia at the onset of the prebreathe procedure. Alveolar 
nitrogen partial pressure (P ly = 11.0 psia) is not used as the initial equilibrium tissue nitrogen par- 
tial pressure to avoid the added complexity of employing the Alveolar Oxygen Equation to estimate 
alveolar nitrogen partial pressure during intermediate altitude exposures where breathing gas 
composition may have been modified. The final nitrogen partial pressure (P t ) in any theoretical 
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tissue type is calculated with equation 2. Ambient nitrogen partial pressure is used to achieve 
operational simplicity in the predicitive models. 

In this example the 360-minute half-time theoretical tissue type is selected. The initial equilibrium 
tissue nitrogen partial pressure (P ) isll.6psia. The prebreathe medium (pure oxygen) contains no 
nitrogen so the pressure difference (P -P 0 ) for the removal of nitrogen is maximal. How rapidly the 
nitrogen is removed from the tissue by this pressure difference from the 360-minute half-time theo- 
retical tissue type is determined by k, the tissue rate constant in equation 2. 

After 4.0 hours of breathing pure oxygen, a final tissue pN 2 (P t ) for the 360-minute half-time tissue 
type is 7.3 psia. In theory, this individual can now decompress from 14.7 psia to 7.3 psia (18,000 feet) 
without supersaturating this theoretical tissue type. If, however, this individual is decompressed to 
4.3 psia the 360-minute TR becomes 


TR = 


7.3 psia 

4.3 psia 


1.70 


If a large group was exposed to 4.3 psia after a 4.0-hour oxygen prebreathe, the true incidence of DCS 
and VGE would eventually be expressed. We can now pair a theoretically derived TR based on a spe- 
cific prebreathe procedure with an experimentally determined percentage of VGE or DCS to obtain a 
data point. 


A different prebreathe procedure, analyzed in the same way, using the 360-minute theoretical half- 
time tissue type, will produce a different TR. This TR matched with VGE or DCS incidence obtained 
from testing the new decompression sickness prevention procedure provides another data point. If 
enough procedures are tested, providing a range of TRs for a specific theoretical half-time tissue, an 
equation can be fitted to the data points thus forming the basis for predicting the outcome of any pre- 
breathe procedure where a TR is calculated. In this way, the TR serves as an empirical value associ- 
ated with some DCS and VGE incidence rather than an absolute value that must not be exceeded. 


DEVELOPMENT OF THE % DCS AND % VGE PREDICTIVE EQUATIONS 


Basic Approach 

To develop the predictive equations it is necessary to determine the theoretical half-time tissue type 
that provides the most precise estimator of subsequent DCS and/or VGE incidence. This is accom- 
plished by accepting the theoretical half-time tissue type that provides the greatest correlation coef- 
ficient after a linear regression is fitted through the DCS or VGE incidence data paired to the calcu- 
lated TRs from a selected theoretical half-time tissue type. A FORTRAN program* that calculates 


*B. F. Edwards, Technology Inc., Houston, Texas, 1986 
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the TRs for a spectrum of theoretical half-time tissue types is available to fascilitate this process. 

It pairs the TR data with the actual DCS and VGE incidence produced under specific prebreathe 
procedures, and performs a linear regression and correlation coefficient analysis of the data. The 26 
prebreathe procedures evaluated by this analysis are outlined in Appendix A. The DCS and VGE 
incidence results from testing the prebreathe procedures are in table 1. 


When a column of TRs from any theoretical half-time tissue type is plotted against DCS or VGE inci- 
dence, the resulting scattergram could indicate a nonlinear relationship between the TRs and the 
% VGE or % DCS data. To develop a mathematical model relating the incidence of VGE and DCS 
to TR two questions had to be answered: (1) What is the best tissue half-time for calculating TRs?, 
and (2) What is the best functional form for the mathematical model? Analysis in the next two 
sections of this report address the questions. 


Fitting Linear Equations to % DCS and % VGE Data 

The first step in developing the predictive equations was to fit a linear function 6 to DCS and VGE 
incidence data paired to calculated TRs from selected theoretical half-time tissue types. A measure 
of the "goodness of fit" of the function to the data was provided by r, the correlation coefficient for the 
regression. An analysis of the correlation data showed that the linear function (straight line) through 
the data provided an acceptable mathematical fit to the data. Figure 1 shows the correlation coeffici- 
ents (r) derived from a linear analysis plotted against the spectrum of theoretical half-time tissue 
types covering from 60 to 960 minutes. 

The point on the curve that indicates the maximum correlation between the % VGE or °h DCS and 
the TR data was used to identify the theoretical half-time tissue type that could be used in predictive 
equations. Table 1 shows the theoretical half-time tissue types, and their TR's paired with actual 
% DCS and °h VGE data. The 360-minute theoretical half-time tissue type TRs show the best linear 
correlation to the % VGE data, while the 480-minute theoretical half-time tissue type TRs show the 
best linear correlation to the % DCS data. 

In the operational use of the models, both DCS and VGE TRs are standardized on the 360-minute 
theoretical half-time tissue type.The predictive linear equations obtained by regressing experimental 
DCS and VGE incidence onto computed TR's, standardized on the 360-minute theoretical half-time 
tissue are 


%VGE = 95.5*Tissue Ratio 

(360') - 

101.1 

r = .88 

(4) 

% DCS = 57.9*Tissue Ratio 

(360') - 

74.1 

r = .78 

(5) 


The predictions from equations 4 and 5 are based on EVA's up to 6.0 hours in duration. Since expo- 
sure time at a final suit pressure is a critical factor in determining DCS or VGE outcome, these equa- 
tions are useful in predicting results from EVA's of 3.0 to 6.0 hours in duration. The section entitled 
"Time to Onset of DCS and VGE During Simulated EVA” addresses DCS and VGE incidence as a 
function of EVA exposure time. 
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TABLE l.-% DCS AND % VGE DATA AND CORRESPONDING TISSUE RATIO 
VALUES COMPUTED AT THE OPTIMUM THEORETICAL TISSUE NITROGEN 
HALF-TIME FOR 34 ALTITUDE DECOMPRESSION EXPERIMENTS 


Prebreathe* 

Procedure 

Number 

%DCS 

480 

Min. 

TR 

%VGE 

360 

Min. 

TR 

Expos. 

Time 

(hrs) 

Expos. 

Press. 

(psia) 

N 

Sex 

Testing 

Inst. 

1 

36 

1.99 

64 

1.80 

3 

4.3 

11 

M 

NASA 

2 

31 

1.96 

89 

1.83 

3 

4.3 

16 

M 

NASA 

3 

41 

1.83 

66 

1.67 

3 

4.3 

12 

M 

NASA 

4 

30 

1.91 

65 

1.70 

4 

4.3 

23 

M 

NASA 

5 

27 

1.89 

45 

1.75 

4 

4.3 

22 

M 

NASA 

6 

21 

1.83 

46 

1.60 

6 

4.3 

28 

M 

NASA 

7 

21 

1.58 

36 

1.45 

6 

4.3 

14 

M 

NASA 

8 

20 

1.82 

57 

1.68 

6 

4.3 

35 

M 

NASA 

9 

16 

1.45 

42 

1.43 

6 

4.3 

12 

M 

NASA 

10 

8 

1.82 

58 

1.68 

3 

4.3 

12 

M 

NASA 

11 

0 

1.31 

17 

1.12 

3 

4.3 

12 

M 

NASA 

12 

0 

1.39 

33 

1.37 

3 

4.3 

12 

M 

NASA 

13 

0 

1.27 

0 

0.95 

3 

4.3 

12 

M 

NASA 

14 

0 

1.33 

33 

1.35 

3 

4.3 

12 

M 

NASA 

15 

0 

1.01 

0 

0.94 

3 

4.3 

12 

M 

NASA 

16 

3 

1.49 

56 

1.49 

6 

7.8 

32 

M 

USAF 

16a 

3 

1.49 

31 

1.49 

6 

7.8 

32 

F 

USAF 

17 

0 

1.40 

45 

1.43 

6 

7.8 

31 

M 

USAF 

17a 

6 

1.40 

35 

1.43 

6 

7.8 

31 

F 

USAF 

18 

0 

1.39 

58 

1.42 

6 

7.8 

31 

M 

USAF 

18a 

10 

1.39 

24 

1.42 

6 

7.8 

29 

F 

USAF 

19 

0 

1.29 

22 

1.29 

6 

9.0 

16 

M 

USAF 

20 

0 

1.16 

13 

1.16 

6 

10.0 

8 

M 

USAF 

21 

0 

1.36 

16 

1.36 

6 

8.5 

9 

M 

USAF 

22 

0 

1.22 

8 

1.22 

6 

9.5 

6 

M 

USAF 

22a 

0 

1.22 

0 

1.22 

6 

9.5 

11 

M 

USAF 

22b 

0 

1.22 

20 

1.22 

6 

9.5 

20 

F 

USAF 

23 

5 

1.58 

32 

1.32 

6 

4.3 

19 

M 

NASA 

23a 

16 

1.58 

26 

1.32 

6 

4.3 

19 

F 

NASA 

24 

0 

1.04 

0 

1.05 

6 

4.3 

8 

M 

NASA 

25 

7 

1.27 

13 

1.24 

6 

6.0 

15 

M 

NASA 

25a 

0 

1.27 

7 

1.24 

6 

6.0 

14 

F 

NASA 

26 

5 

1.40 

25 

1.40 

6 

8.3 

20 

M 

USAF 

26a 

0 

1.40 

27 

1.40 

6 

8.3 

11 

F 

USAF 


N = 607 


♦See Appendix A for description of each Prebreathe Procedure. 
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Figure 1.- Plots of the linear correlation coefficient versus specific theoretical half-time tissue types for % DCS and % VGE data. 


Equations 4 and 5 were derived from a subset of the data from table 1. Not all the data on % DCS 
and % VGE contained in table 1 were used to develop these equations. The reason for excluding 
some data is that prebreathe procedures that provided complete bends protection also produced zero 
incidence of DCS and VGE. TRs derived from an analysis of these prebreathe procedures are also 
very low. As a result, incorporating these zero % DCS and % VGE data will artificially shift the slope 
of the linear equation. Excluding some zero data eliminates shifting the slope of the line to include 
prebreathe conditions that were extremely conservative. If some zero incidence data points are 
excluded, the predictive equations provide a better data fit. 

The excluded zero data consisted of those zero percent DCS and VGE data points that had 360-min- 
ute theoretical half-time TRs that were less than the TR that first produced some DCS or VGE inci- 
dence. Using this method, the % DCS data set included 26 data points and the % VGE data consisted 
of 31 points (see table 1). 

Figures 2, 3 and 4 show the linear plots of % DCS and % VGE against the TRs calculated from the 
360-minute theoretical half-time tissue type. The 95% confidence band for figures 2 and 3 is 
included to quantify the variability in the DCS and VGE data. 

These equations can predict the best estimate of DCS or VGE incidence within specified confidence 
limits. The % DCS linear regression indicates a zero incidence of Type I DCS when a TR of 1.28 is 
achieved. The % VGE linear regression indicates a zero incidence of VGE when a TR of 1 .06 is 
achieved. However; the incidence of DCS and VGE in figures 2 and 3 indicate that the incidence 
of DCS and VGE at these TRs can be greater than zero. To more accurately describe the relationship 
between TR and % DCS and % VGE when the TR is very low ( < 1 .2) or very high ( > 1.8) the Hill 
equation is employed. 


Fitting the Hill Equation to % DCS and % VGE Data 

Since the relationship of both the incidence of DCS and VGE to TR are bounded at 0% and 100%, a 
model having sigmoidal characteristics is suggested. The two parameter Hill equation provides a 
model for such a sigmoidal relationship. The Hill equation, a probabilistic dose-response function, is 
frequently used in biology and pharmacology to quantify the dose-response behavior of subjects to 
some stress, medication, or treatment. The Hill equation has the form 

P = D n /(D n +D^) (6) 

where D is the dose, P is the probability of a positive response in a dichotomous (yes/no) outcome set, 
and D 50 is the dose at which P equals 0.5 (50%). The exponent n and D 50 are parameters to be esti- 
mated from the experimental data. Equation 6 is sigmoidal in shape, having a value of zero for zero 
dose, a value of 0.5 at and, approaches unity as the certainty of a positive response increases 
with increasing dosage. The exponent n, the order of the Hill equation, controls the slope of the 
central portion of the curve as it passes through the 50% point. 
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Figure 2.- Linear plot of DCS incidence using 360-minute theoretical half-time tissue and its Tissue Ratios 

with the 95% confidence interval. 
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360 Minute Tissue Ratio 

% Symptoms Resulting in Test Termination plus Delayed Symptorr.s= 10.5* Tissue Ratio (360')- 13.1 r 

Figure 4 - Linear plot of DCS symptoms that resulted in test termination or symptoms that reoccurred after the test. 


For probabilistic modeling of DCS, VGE, or the more serious DCS symptoms where TR quantifies the 
dose, the Hill equation is modified to reflect an incidence probability of zero at a TR of 0.78. This is 
required since the original definition of TR developed in equation 3 calculates a TR of 0.78 under 
equilibrium conditions at sea level. A TR of 0.78 becomes the starting-point or base line from which 
to construct the Hill equation. Therefore, the more appropriate form of the Hill equation in this 
application is 

TR- .78 j "/ 

The parameter n and TR 50 , and the correlation coefficient (r) of the fit are estimated using the ex- 
perimental data set {(TR, P)} found in table 1, and a proprietary non-linear estimation program 14 . 
Unlike the zero exclusion procedure used in the linear regressions, the Hill equation utilizes all data 
from table 1. 

When fitted to the experimental data the Hill equations for predicting % VGE, % DCS, and % DCS 
Refractory Symptoms are 




and 


/ \ 3-08 

\( 


>. 3.08 




% VGE = (TR- .78 j 

/ 

(TR- 

.78 

) + 0.47 

*100 

r = .91 

(8) 

/ \ 4-24 | 

7 

\ 

, 4.24 




%DCS = (JTR- .78 j 

/ 

l( TR - 

.78 

> 

J +2.16 

*100 

r = .84 

(9) 


% Refractory Symptoms 


/ \ 2-50 

[/ \ 2.50 1 

(TR— .78 j / 

{ TR - .78 j +17.61 


*100 r = .53 


where TR is the Tissue Ratio for the 360-minute theoretical half-time tissue, and % Refractory 
Symptoms (Grade 3) are those DCS symptoms that resulted in test termination or symptoms 
that reoccurred after the test. 


GO) 


Figures 5, 6 and 7 display the Hill equations with the experimental data points and the 95% confi- 
dence interval. Figure 8 shows a composite of the three curves. Figures 9 and 10 show each curve 
without the 95% confidence interval. These curves, in addition to the curve in figure 7, are currently 
the accepted predictors of DCS and VGE incidence based on all available data that can be applied to 
current Shuttle EVA operations. 


14 
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Figure 5.- Hill plot of DCS incidence using 360-minute theoretical half-time tissue and its Tissue Ratios with the 95% confidence interval. 
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Figure 6.- Hill plot of VGE incidence using 360-minute theoretical half-time tissue and its Tissue Ratios with the 95% confidence interval 
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Figure 7.- Hill plot of DCS symptoms that resulted in test termination or symptoms that reoccurred after the test. 
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Figure 8.- Hill plot of DCS and VGE incidence using 360-minute theoretical half-time tissue and its Tissue Ratios 
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Figure 9.- Hill 


Hill 's Equation: 

% VGE = (TR - .78) 3 0Q / [ (JR - ,78) 3 ' 08 + .47] * 100 
r = 0.91 n = 34 

Based on 607 human decompressions with LVA 
work simulation 
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Figure 10.- Hill plot of VGE incidence using 360-minute theoretical half-time tissue and its Tissue Ratios. 



Confidence Intervals for the Predictive Linear and Hill Equations 


Associated with each prediction of % DCS or % VGE is a degree of variability related to the scatter 
of real data about the linear regression line. A measure of this scatter is the error variance for the 

regression. From Ostle 5 , the error variance j is given by 

s '=l( Y - ? ) 2 '("- 2 ) 

where Y is the measured DCS or VGE percentage value, Y is the predicted value (a function of the 
TR and provided by the predictive model), and n is the number of data points in the real data set. 


The estimated variance 


^S^j for a predicted group outcome for a TR value x is 


( - V 2 

t \ 2 i 

1 + (1/n) + 1 x — x J 

l*i-V 1 


where the are the TR data set, and x is the mean TR for the data set. 

Finally the confidence interval (Cl) for the DCS or VGE incidence estimate at TR x is given by 

CI= Y ± t (1 _ a/2) ( n _ 2 ) * S f 

where t is the specified t-statistic, and (1-a) is the desired level of confidence, i.e., 95%. 

For an example, the linear equation 4 predicts a 42.0% incidence of VGE for a TR of 1.50 calculated 
from the 360 minute theoretical half-time tissue (fig. 3). The 95% confidence interval calculated for 
this prediction is 42.0% ± 20.0%. Since the true VGE incidence at a TR of 1.50 is a fixed value rather 
than a random variable, an interpretation of the confidence interval is that one can be 95% certain 
that it will contain the true VGE incidence. Alternately, 95% of the confidence intervals computed 
from similar data samples will contain the true VGE incidence. 

Equation 5 predicts a 13.0% incidence of DCS with a TR of 1.50 calculated from the 360-minute 
theoretical half-time tissue type. The 95% confidence interval for this prediction produces an upper 
limit of 30.0% and a lower limit of 0.0% incidence within a tested group. 

This same procedure can be applied to any prediction of % DCS or % VGE when using a spectrum 
of TRs. Figures 2 and 3 show a 95% confidence band about the linear regression line of % DCS and 
% VGE against their respective TRs. Thus, given a specific TR, a range of predicted % VGE or 
%DCS incidence can be determined. 

The Hill equations (eq. 8 and 9) are analyzed in the same fashion. Figures 5 and 6 show the 95% con- 
fidence interval about the sigmoidal functions. In figure 5, a TR of 1 .50 would predict a 10% DCS 
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incidence. The 95% confidence interval calculated for this prediction is 10% ± 14%. In figure 6, the 
same TR predicts 44% VGE with a 95% confidence interval of 44% ± 20%. 

Since the data pairs that comprise all the figures were from group results, it must be stressed that the 
95% confidence interval applies to the percent incidence you would expect if a group of individuals, 
randomly chosen from a general population, were tested. For a group, the a priori expected incidence 
and the a posteriori incidence should be close, i.e., the a posteriori incidence should be within the con- 
fidence interval. For an individual, an a priori level of risk can be predicted, however; the a posteriori 
incidence will be 0% or 100%. 


Accuracy and Precision of Predictive Linear and Hill Equations 

Accuracy assessment of equations 4,5,8, and 9 for predicting the incidence of DCS and VGE is accom- 
plished using a cross validation procedure . The DCS and VGE data set is randomly divided into two 
equal subsets. One subset is designated the training set, and is used to generate a new linear and Hill 
predictive equations. The remaining data for this cross validation test is the test set . 

Substituting the test set TR values into the training set predictive equations provides a set of % VGE 
and % DCS values and residuals which are statistically compared to the observed percentage values 
and residuals of the training set. Residuals are the differences between observed values and the 
regression equation predicted values. 

Table 2 contains the training set and test set statistics for both the linear and Hill equation data, i.e., 
the average (d) and standard deviation (s d ) of the unsigned residuals, the mean error of the fit (e). 
the 95% confidence upper bound for the test set deviation, and the critical and computed F-statistics 
for testing the hypothesis that the mean square residuals of the training set and the average squared 
residuals of the test set are equal. 

From table 2, F-statistic comparisons provide no basis for rejecting the null hypothesis that the preci- 
sion with which the training set predictive equations represent the test set is significantly different 
from that with which the predictive equation represents the training set. However, the mean devia- 
tions of the test sets are larger than those of the training sets except for the % DCS Hill equation. In 
regard to the average deviation upper bound, loosely interpreted, one can be 95% confident that the 
incidence of VGE or DCS will be no more than approximately 20% greater than that predicted by 
equations 4,5,8, and 9. This is the degree of precision possible for any % VGE or % DCS prediction. 
Finally, in regard to accuracy, the average error (e) for each test fit is small, 1.0% or less, indicating 
the central location of the predictive equations with respect to both training and test sets. The excep- 
tion is the % VGE Hill equation. 


TIME TO ONSET OF DCS AND VGE DURING SIMULATED EVA 

In the preceding analysis, consideration was given to the occurrence of DCS or VGE without regard 
to the time of onset or detection after the beginning of the simulated EVA. In this section the com- 
bined Air Force 3,4,17,18 and NASA data 7,8,9,10,11,15 ’ 16 are analyzed to provide graphs for predicting 
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TABLE 2.- CROSS VALIDATION TRAINING AND TEST STATISTICS 


Linear Equations 


VGE Data 

N 

*.95(n — 1) 

d 

s a 

df 

.95 

e 

F .975(dfl,df2) F 1 

r*) /(*). 

\ d/max y d^min 

Test Set 

15 

1.761 

8.13 

6.17 

15 

19% 

1.06% 

2.86 

1.39 

Trng. Set 

16 


8.06 

5.58 

15 


0.38% 



DCS Data 










Test Set 

13 

1.782 

8.62 

5.06 

13 

18% 

-0.92% 

3.31 

1.26 

Trng. Set 

13 


4.62 

4.50 

12 


0.08% 




Hill Equations 


VGE Data 

N 

*.95(11- 1) d 

s d 

df 

.95 

7 

F F 1 

.975<dfl t df2) 1 

( s a) '( 

\ dymax \ 

i) 

d ymin 

Test Set 

17 

10.94 

6.55 

17 


4.47% 






1.746 



22% 


2.70 

1.05 


Trng. Set 

17 

4.88 

6.24 

16 


0.06% 




DCS Data 










Test Set 

17 

2.59 

3.20 

17 


-0.47% 






1.746 



8% 


2.70 

2.68 


Trng. Set 

17 

6.41 

8.58 

16 


-0.29% 
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the incidence of DCS and VGE as a function of the time into the EVA at a specified TR. The data in 
Table 1 consists of 1 1 1 tests for 3 hours of simulated EVA exposure, 45 tests for 4 hours of exposure, 
and 451 tests for 6 hours of EVA exposure. 

Even though the expected DCS or VGE incidence may be quite high at certain TRs, the risk can be 
considerably lessened if the EVA time is shortened. Figures 11 and 12 present the cumulative inci- 
dence of DCS and VGE as linear functions of TR at various simulated EVA times up to 360 minutes. 

The time lines in figures 1 1 and 12 are structured as follows 

• Equations 4 and 5 are solved for the TR intercepts to obtain 
% VGE = 0 at TR = 1.057, and 

% DCS = 0 atTR = 1.279 

• Assuming that all cumulative DCS or VGE time lines pass through the same respective TR 
intercept, the cumulative DCS and VGE incidence and TR data at the specified time are fitted 
to respective linear functions constrained to pass through the corresponding TR intercept. 
Figures 1 1 and 12 constitute families of cumulative percent incidence versus TR lines that 
have time as a variable parameter. 

Figures 13 and 14 are families of lines for DCS and VGE incidence versus time with TR as a variable 
parameter. These TR lines are obtained by solving each of the family of DCS and VGE time lines for 
the incidence at a given TR and drawing a smooth curve through these calculated data points. Both 
figures show the expected incidence given a particular 360 minute theoretical half-time TR and a 
planned EVA duration. They express the increasing risk of DCS and VGE given a particular TR 
as EVA duration is increased. 

Table 3 contains the calculated linear equations that can be applied to predict % DCS or % VGE for 
a specified EVA duration or if an emergency EVA of limited duration is needed. 


DCS and VGE Latency Times 


Table 4 contains data on the average time to report Type I DCS symptoms or to detect VGE after a 
particular prebreathe procedure. Figure 15 displays the DCS latency times against the calculated TR 
from the 360 minute theoretical half-time tissue type. A nonlinear fitting function was applied to the 
data. Figure 16 displays the VGE latency times plotted against the calculated TR. In both cases, the 
data indicates that lower TRs will delay the initial appearance of VGE and complaints of Type I DCS. 
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Figure 11.- Selected linear equations to describe DCS incidence as a function of simulated EVA exposure time. 
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Figure 12.- Selected linear equations to describe VGE incidence as a function of simulated EVA exposure time. 



o 


■*-> 

CC 


(D 

3 

00 



t 

( 

( 


EVA Exposure Time (Hours) 


2.0 Tissue Ratio 



aouapLoui 3 9 A % 


28 


Figure 14.- 


TABLE 3.- DCS AND VGE INCIDENCE AS A FUNCTION OF EVA DURATION 


EVA Duration 
(minutes) 



Linear Predictive 
Equation 


Correlation 

Coefficient 

N 

30 

%DCS 

— 

5.9 * 

Tissue Ratio (360') - 

7.1 

.10 

6 


%VGE 

= 

26.7 * 

Tissue Ratio (360') - 

34.0 

.63 

31 

60 

%DCS 

= 

21.9 * 

Tissue Ratio (360') - 

29.9 

.74 

23 


%VGE 

= 

58.5 * 

Tissue Ratio (360') - 

72.5 

.77 

31 

90 

%DCS 

= 

49.1 * 

Tissue Ratio (360') - 

68.0 

.82 

23 


%VGE 


76.5 * 

Tissue Ratio (360') - 

91.2 

.74 

31 

120 

%DCS 


54.6 * 

Tissue Ratio (360') - 

74.6 

.80 

23 


%VGE 

= 

88.8 * 

Tissue Ratio (360') - 

104.0 

.79 

31 

150 

%DCS 

= 

47.0 * 

Tissue Ratio (360') - 

62.4 

.77 

26 


%VGE 

= 

103.0 * 

Tissue Ratio (360') - 

120.0 

.87 

31 

180 

%DCS 

zz 

55.7 * 

Tissue Ratio (360') - 

73.7 

.79 

26 


%VGE 

= 

104.0 * 

Tissue Ratio (360') - 

119.0 

.87 

31 

210 

%DCS 


41.7 * 

Tissue Ratio (360') - 

54.3 

.71 

23 


%VGE 

= 

94.2 * 

Tissue Ratio (360') - 

105.0 

.83 

28 

240 

%DCS 

— 

43.0 * 

Tissue Ratio (360') - 

54.6 

.66 

23 


%VGE 


93.5 * 

Tissue Ratio (360') - 

102.0 

.82 

28 

270 

%DCS 

= 

29.5 * 

Tissue Ratio (360') - 

34.8 

.43 

18 


%VGE 

= 

109.0 * 

Tissue Ratio (360') - 

124.0 

.80 

22 

300 

%DCS 

— 

32.2 * 

Tissue Ratio (360') - 

38.4 

.45 

18 


%VGE 

= 

109.0 * 

Tissue Ratio (360') - 

122.0 

.82 

22 

330 

%DCS 


36.0 * 

Tissue Ratio (360') - 

43.6 

.49 

18 


%VGE 

= 

107.0 * 

Tissue Ratio (360') - 

119.0 

.83 

22 

360 

%DCS 

= 

34.0 * 

Tissue Ratio (360') - 

40.3 

.48 

18 


%VGE 


104.0 * 

Tissue Ratio (360') - 

114.0 

.82 

22 
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TABLE 4.- DCS AND VGE LATENCY TIMES 


Prebreathe Average VGE Average DCS 

Procedure 360 Final Exposure Latency Time Latency Time 

Number TR Pressure (psia) (minutes) ± (SD) N (minutes) + (SD) N 


1 

1.80 

4.3 

54 

(23) 

7 

61 

(15) 

4 

2 

1.83 

4.3 

47 

(28) 

14 

82 

(52) 

5 

3 

1.67 

4.3 

48 

(31) 

7 

109 

(37) 

4 

4 

1.70 

4.3 

120 

(62) 

15 

105 

(60) 

7 

5 

1.75 

4.3 

43 

(40) 

10 

134 

(68) 

6 

6 

1.60 

4.3 

60 

(41) 

13 

132 

(84) 

6 

7 

1.45 

4.3 

140 

(58) 

5 

165 

(81) 

7 

8 

1.68 

4.3 

124 

(70) 

20 

173 

(85) 

3 

9 

1.43 

4.3 

121 

(27) 

5 

228 

(12) 

2 

10 

1.68 

4.3 

105 

(43) 

7 

24 

(-) 

1* 

11 

1.12 

4.3 

45 

(17) 

2 * 




12 

1.37 

4.3 

88 

(45) 

4 




14 

1.35 

4.3 

66 

(58) 

4 




16/1 6a 

1.49 

7.8 

108 

(64) 

28 




17/1 7a 

1.43 

7.8 

122 

(63) 

25 




18/1 8a 

1.42 

7.8 

101 

(62) 

24 




19 

1.29 

9.0 

96 

(38) 

7 




20 

1.16 

10.0 

242 

(6) 

2 




21 

1.36 

8.5 

195 

(81) 

3 




22/22a 

1.22 

9.5 

293 

(44) 

4 




23/23a 

1.32 

4.3 

211 

(99) 

11 

142 

(68) 

4 

25/25a 

1.24 

6.0 

150 

(59) 

3 

150 

(-) 

1* 

26/26a 

1.40 

8.3 

118 

(98) 

8 





♦Data omitted from mathematical analysis shown in figures 15 and 16. 
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Figure 15.- Latency time for initial DCS symptoms versus 360-minute theoretical half-time Tissue Ratios 





300 



(* uilu) auji_L Adu3}D“| 39 A '& A V 


32 


Figure 16.- Latency time for initial VGE detection versus 360-minute theoretical half-time Tissue Ratios. 



SUMMARY OF DCS AND VGE PREDICTIVE EQUATIONS 


Tissue Nitrogen Partial Pressure Relaxation Equation 


P .= I> „ + ( P .- P o)( 1 -‘>' k, ) 

where 

P = the nitrogen partial pressure in the tissue after exposure for t minutes, 

P Q = initial tissue nitrogen partial pressure, 

P g = ambient nitrogen partial pressure in breathing medium, 
e = base of natural logarithm, 
t = exposure time in minutes, and 
k = tissue nitrogen partial pressure rate constant. 

The tissue rate constant (k) is related to the tissue nitrogen half-time (t 1/2 ) through the equation, 


where 


Tissue Ratio 


k = (ln2) / tj_ 2 = -693/t^ 


= tissue nitrogen partial pressure half-time 


Final Tissue pN (abs) 

TR = 

Ambient Pressure (abs) 
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Linear DCS and VGE Predictive Equations 


%VGE = 95.5* Tissue Ratio (360‘)- 101.1 r=.88 

% VGE = 57.9 * Tissue Ratio (360* ) - 74. 1 r = .78 
Hill DCS and VGE Predictive Equations 


/ \ 3.08 r / \3.UB 

% VGE = f TR- .78 j /MTR-.78J +0.47 

%DCS = (tR— .78 j / | ^TR — .78^ +2.16 


*100 r = .91 


*100 r = .84 


/ \ 2.50 r / \2.50 

% Refractory Symptoms = (TR— .78 j / ^TR — .78 j +17.61 


100 


where TR is the Tissue Ratio of the 360 minute theoretical half-time tissue. 



USER SECTION 


Predicting the Incidence of DCS and VGE 


In order to properly calculate TRs for subsequent use in equations 8,9, and 10 the user must follow the 
exact procedures for calculating TRs that were used in the development of the predictive equations. 

A simple example of the use of equations 2 and 3 was given previously. A more complex example is 
presented to aquaint the user with three important rules that must be understood when calculating 
the 360 minute theoretical TR. A second example on the proper use of the models during a simulated 
emergency EVA is presented to demonstrate how to design a prebreathe procedure that takes 
advantage of a short EVA duration. 


Example 1 

Problem Calculate the risk of Type I DCS, VGE, and DCS symptoms that could result in early 

EVA termination given the following prebreathe procedure: 

Event 1 Two crewmen prebreathe 100% oxygen for 60 minutes while at 14.7 psia in the Shuttle. 

Event 2 The Shuttle atmospheric pressure is reduced to 10.2 psia within 10 minutes, while the 
EVA crewmen are still breathing 100% oxygen. Breathing gas composition at 10.2 psia 
is adjusted to 27% oxygen, 73% nitrogen. The entire crew spends 24 hours in this new 
pressure/gas environment. 


Event 3 A 40 minute, 100% oxygen prebreathe is performed in the suit while still at 10.2 psia 
cabin pressure. This is followed by a 30 minute decompression to 4.3 psia. 

Problem Calculate the resulting 360 minute theoretical half-time TR using equations 2 and 3 
for the first EVA. Use the calculated TR in conjunction with equations 8,9, and 10 to 
calculate the risk of DCS and VGE for the first EVA. 


Event 4 4.0 hours of EVA are performed while breathing 100% oxygen at 4.3 psia. The crew 

then returns to the Shuttle for a 2.0 hour lunch. Shuttle atmosphere is at 10.2 psia 
with 27% oxygen, 73% nitrogen. 

Event 5 Once again, a 40-minute, 100% oxygen prebreathe is performed in the suit prior to a 
30-minute decompression to 4.3 psia for a 3.0 hour EVA. 

Problem Calculate the resulting 360 minute theoretical half-time TR using equations 2 and 3 
for the second EVA. Use the calculated TR in conjunction with equations 8,9, and 10 
to calculate the risk of DCS and VGE on the second EVA. 
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Example 1 Calculations 
Initial Comments 

(1) The 360 minute theoretical half-time tissue is to be used in all applications of equation 2. 

(2) Time is expressed in minutes. 

(3) Pressure is expressed in psia. 

Rule 1 Always use the ambient nitrogen partial pressure (P a ) in the breathing mixture as the 
driving pressure in equation 2. Do not use alveolar nitrogen partial pressure. Initial 
tissue nitrogen partial pressure (P 0 ) is equal to 11.6 psia if the crewman has been 
exposed to sea level pressure and gas composition for over 24 hours. 

Step 1 Two crewmen prebreathe 100% oxygen for 60 minutes while at 14.7 psia in the Shuttle. 

In equation 2, P Q = 11.6 psia, while P a = 0.0 psia since 100% oxygen contains 
no nitrogen. Equation 2 calculates a final tissue nitrogen partial pressure (P t ) 
of 10.33 psia for the 360-minute theoretical half-time tissue after 60 minutes of 
prebreathing. 

Rule 2 Any ambient pressure transition (decompression or recompression) that exceeds 
10 minutes in duration must be accounted for in the prebreathe procedure. Since 
crewmen are usually breathing 100% oxygen during these transitions, the time from 
one ambient pressure to another must be included as part of the 100% oxygen prebreathe 
period if it exceeds 10 minutes. 

Step 2 The Shuttle atmospheric pressure is reduced to 10.2 psia within 10 minutes while the 
EVA crewmen are still breathing 100% oxygen. Breathing gas composition is adjusted 
to 27% oxygen and 73% nitrogen. The entire crew spends 24 hours in this new 
pressure/gas environment. 

Since the decompression time from 14.7 psia to 10.2 psia was less than 10 minutes, the 
additional oxygen prebreathe time is not added to the initial 60 minute prebreathe. 

The final tissue nitrogen pressure is still 10.33 psia upon reaching 10.2 psia. 

Equation 2 is used again to establish the final tissue nitrogen partial pressure after 
24 hours (1440 minutes) of exposure to 73% nitrogen. In this case 

• P Q = P t from the previous application of equation 2: 10.33 psia. 

• P a = 73% nitrogen * 10.2 psia = 7.44 psia nitrogen. 

• P t after 1440 minutes in this environment = 7.62 psia nitrogen. 
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Step 3 A 40-minute, 100% oxygen prebreathe is performed in the suit followed by a 30-minute 
decompression to 4.3 psia. 

Rule 2 requires that the 30 minutes of oxygen prebreathing during the decompression 
from 10.2 psia to 4.3 psia be included in the oxygen prebreathe period. Total oxygen 
prebreathe = 70 minutes. 

Equation 2 is used again to establish the final tissue nitrogen partial pressure after the 
70-minute exposure to 0.0 psia nitrogen. In this case 

• P q = P t from the previous application of equation 2: 7.62 psia. 

• P a = 0% nitrogen* 10.2 psia = 0.0 psia nitrogen. 

• P t after 70 minutes in this environment = 6.66 psia nitrogen. 

For the first EVA, the TR is calculated by using equation 3 
TR = P t /Ambient Pressure (suit pressure) 

= 6.66 psia/4.3 psia = 1.55 

Substituting this TR into equations 8,9, and 10 results in a 49% risk of VGE, 13% risk of 
Type I DCS, and a 3% risk of DCS symptoms that may force a crewman to end his EVA. 
These predictions are based on EVA durations up to 6.0 hours. 

Rule 3 During multiple exposures to low ambient pressures separated by hours or days at 

higher nitrogen pressures, one must account for the nitrogen removed during the 100% 
oxygen prebreathe while at low pressures. In other words, the time during EVA must be 
used as part of the prebreathe period if additional EVAs are to be planned within 24 
hours. 

Step 4 4.0 hours of EVA are performed while breathing 100% oxygen at 4.3 psia, and the 

crew returns to the Shuttle for a 2.0 hour lunch. 

Equation 2 is used again to establish the final tissue nitrogen partial pressure after 
4.0 hours (240 minutes) of exposure to 0.0 psia nitrogen during the EVA. In this case 

• P = P from the previous application of equation 2: 6.66 psia nitrogen 

O l 

• P a = 0.0 psia nitrogen 

• P after 240 minutes in this environment = 4.20 psia nitrogen. 

The crew returns to the Shuttle environment where P a = 7.44 psia for 2.0 hours 
(120 minutes). P 0 = Pt = 4.20 psia nitrogen and equation 2 calculates that P t 
after 2.0 hours in the Shuttle is 4.86 psia nitrogen. 
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Step 5 Once again, a 40-minute oxygen prebreathe is performed in the suit prior to a 30 minute 
decompression to 4.3 psia for a 3.0-hour EVA. 

Rule 2 also applies in this situation. Seventy minutes of 100% oxygen prebreathe must 
be used in the calculation. 

Equation 2 is used again to establish the final tissue nitrogen partial pressure after 
70 minutes of exposure to 0.0 psia nitrogen. In this case 

• P = P from the previous application of equation 2: 4.86 psia nitrogen 

O l 

• P = 0.0 psia nitrogen 

• P t after 70 minutes in this environment = 4.25 psia nitrogen. 

For the second EVA, TR is calculated using equation 3 

TR = Pt /Ambient Pressure (suit pressure) 

= 4.25 psia/4.3 psia = 0.98 

Substituting this TR into equations 8,9, and 10 results in a 2% risk of VGE, 0% risk 
of Type I DCS, and a 0% risk of DCS symptoms that could terminate an EVA. 


Example 2 

Problem Shuttle crew members are at 14.7 psia breathing 2 1% oxygen and 79% nitrogen 

atmosphere. A life-threatening failure in the cargo bay occurs that requires an EVA 
for a repair mission. EVA duration should not exceed 2.0 hours for this repair task. 
Calculate the minimum time of 100% oxygen prebreathing that would allow a safe 2.0- 
hour EVA. In this emergency, a safe EVA is defined as one that does not expose EVA 
crewman to a Type I DCS risk greater than 25%. 

From equation 9, a TR of 1 .7 would predict a 25% risk of Type I DCS. Equations 4, 5, 8, 9, 
and 10 were derived from the cumulative incidence of DCS and VGE from simulated 
EVAs up to 6.0 hours. The incidence of DCS and VGE they predict reflect the influence 
of long duration EVAs. In this simulated emergency situation where EVA duration is 
planned for 2.0 hours, a higher TR can theoretically be used that provides the same DCS 
risk if the EVA exposure is limited. Figure 13 graphically shows that a TR of 1.9 could 
be tolerated for 2.0 hours and still predict a 25% risk of Type I DCS. In addition, Figure 
15 shows that the average time to detect DCS with a 1.9 TR is 110 minutes. This would 
indicate a 2.0 hour EVA would present an acceptable risk for this emergency. Equation 
2 can now be evaluated such that 100% oxygen prebreathe time becomes the variable to 
calculate, while the final TR becomes the fixed value. In this situation, the least amount 
of oxygen prebreathing prior to decompression to 4,3 psia is 182 minutes . 
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If the emergency EVA condition warrants less prebreathe time, then the risk of Type I 
DCS (higher TR) could be increased, or the planned EVA duration could be shortened. 
This flexibility allows the user options in the event of emergency EVA conditions. 


Assumptions in and Limitations of The Models 


Scientists and engineers attempt to maximize and utilize man's ability to adapt to low pressure en- 
vironments. Space suits operating at a pressure lower than cabin pressure are desirable for many 
reasons and have been successfully used throughout the U.S. space program. The need to properly 
manage denitrogenation-decompression procedures prior to a reduction in ambient pressure is crit- 
ical in the one atmosphere Shuttle and Space Station environments. The models developed in this 
report are useful to that end, but their limitations are discussed. 


Any empirical model is limited in its ability to predict. It is limited to the data used to derive the 
model. Prebreathe conditions, possibly dictated by emergency conditions, that deviate dramatically 
from those tested may produce outcomes not predicted by the models. Denitrogen-decompression 
procedures that fall outside the range of the tests used to define these models or which deviate in any 
of the parameters that are known or suspected to effect the incidence of DCS may produce outcomes 
not predicted by the models. In such cases it is critical to perform verification testing prior to 
accepting the results from the models. 


One of the environmental parameters that is different in space flight than in the altitude chamber 
test programs is the lack of gravity, and the physiological adaptations it induces. It is possible that 
microgravity could have some influence on the incidence of DCS or VGE. At the present time there 
is no positive or negative evidence of a zero-g influence on altitude decompression sickness. Experi- 
ments being planned for Space Station may provide evidence of a microgravity effect on DCS or 
VGE incidence during EVA. 
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APPENDIX A 

OUTLINE OF 26 SEPARATE PREBREATHE PROCEDURES 
USED IN 607 MANNED TESTS 
CONDUCTED AT THE JOHNSON SPACE CENTER 

AND 

BROOKS AIR FORCE BASE (1982 - 1986) 


No. Prebreathe Procedure 

1. 3.5 hours oxygen prebreathe at 14.7 psia prior to 3.0 hour exposure to 4.3 psia. 
Decompression was rapid. Exercise stressed lower body. N = 1 1 

2. 12.0 hours at 10.2 psia plus a 40 minute oxygen prebreathe prior to a 3.0 hour exposure 
to 4.3 psia. Decompression was rapid. Exercise stressed lower body. Gas composition at 

10.2 psia was 26.5% 0 2 - 73.5% N 2 - N = 16 

3. 12.0 hours at 10.2 plus a 90 minute oxygen prebreathe prior to a 3.0 hour exposure to 
4.3 psia. Decompression was rapid. Exercise stressed lower body. Gas composition at 

10.2 psia was 26.5% 0 2 - 73.5% N 2 . N = 12 

4. 3.5 hours oxygen prebreathe at 14.7 prior to a 4.0 hour exposure to 4.3 psia. Decompression 
was gradual and allowed 30 minutes of additional oxygen prebreathe prior to reaching 

4.3 psia. Exercise stressed upper body. N = 23 

5. 12.0 hours at 10.2 psia plus a 40 minute oxygen prebreathe prior to a 4.0 hour exposure 
to 4.3 psia. Decompression was gradual and allowed 25 minutes of additional oxygen pre- 
breathe prior to reaching 4.3 psia. Exercise stressed upper body. Gas composition at 

10.2 psia was 26.5% 0 2 - 73.5% N 2 . N = 22 

6. 4.0 hours oxygen prebreathing at 14.7 psia prior to a 6.0 hour exposure to 4.3 psia. Decom- 
pression was gradual and allowed 30 minutes of additional oxygen prebreathe prior to 
reaching 4.3 psia. Exercise stressed upper body. N = 28 

7. Same procedure as #6 except crew returned to 14.7 psia for 17.0 hours. Second EVA began 
after 4.0 hours oxygen prebreathe at 14.7 psia prior to second 6.0 hour exposure to 4.3 psia. 
Decompression was gradual and allowed 30 minutes of additional oxygen prebreathe prior 
to reaching 4.3 psia. Exercise stressed upper body. N = 14 

8. 60 minutes oxygen prebreathe at 14.7 psia followed by 12.0 hours at 10.2 psia plus an 
additional 40 minute oxygen prebreathe prior to a 6.0 hour exposure to 4.3 psia. Decom- 
pression was gradual and allowed 25 minutes of additional oxygen prebreathe prior to 
reaching 4.3 psia. Exercise stressed upper body. Gas composition at 10.2 psia was 26.5% 
0 2 - 73.5% N 2 . N = 35 
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9. Same procedure as #8 except crew returned to 10.2 psia for 17.0 hours. Second EVA began 
after 40 minute oxygen prebreathe at 10.2 prior to a second 6.0 hour exposure to 4.3 psia. 
Decompression was gradual and allowed 25 minutes of additional oxygen prebreathe prior 
to reaching 4.3 psia. Exercise stressed upper body. N = 12 

10. 60 minutes oxygen prebreathe at 14.7 psia followed by 12.0 hours at 10.2 psia plus an 
additional 40 minute oxygen prebreathe prior to a 3.0 hour exposure to 4.3 psia. This 
was the first of two exposures in the same day. Decompression was gradual and allowed 
25 minutes of additional oxygen prebreathe prior to reaching 4.3 psia. Gas composition 
at 10.2 psia was 26.5% 0 2 - 73.5% N 2 . N = 12 

11. Same procedure as #10. Crew then returned to 10.2 psia for 80 minutes. A 40 minute 
oxygen prebreathe was then performed prior to a second 3.0 hour exposure to 4.3 psia in 
the same day. Decompression was gradual and allowed 25 minutes of additional oxygen 
prebreathe prior to reaching 4.3 psia. N = 12 

12. Same procedure as #10 plus #1 1 except crew returned to 10.2 psia for 14.0 hours. First 
EVA of second day began with a 40 minute oxygen prebreathe prior to a 3.0 hour exposure 
to 4.3 psia. Decompression was gradual and allowed 25 minutes of additional oxygen 
prebreathe prior to reaching 4.3 psia. N = 12 

13. Same procedure as #8 plus #11 plus #12. Crew then returned to 10.2 psia for 80 minutes, 
a 40 minute oxygen prebreathe was then performed prior to a second 3.0 hour exposure of 
the second day EVA tp 4.3 psia. Decompression was gradual and allowed 25 minutes of 
additional oxygen prebreathe prior to reaching 4.3 psia. N = 12 

14. Same procedure as #10 plus #11 plus #12 plus #13 except crew returned to 10.2 psia for 

14.0 hours. First EVA of third day began with a 40 minute oxygen prebreathe prior to a 

3.0 hour exposure to 4.3 psia. Decompression was gradual and allowed 25 minutes of 
additional oxygen prebreathe prior to reaching 4.3 psia. N = 12 

15. Same procedure as #10 plus #1 1 plus #12 plus #13 plus #14. Crew then returned to 
10.2 psia for 80 minutes. A 40 minute oxygen prebreathe was then performed prior to a 
second 3.0 hour exposure of the third day EVA to 4.3 psia. Decompression was gradual 
and allowed 25 minutes of additional oxygen prebreathe prior to reaching 4.3 psia. 

N = 12 

16. Exposure to 7.8 psia for 6.0 hours using 50% O 2 - 50% N 2 mixture without prior oxygen 
prebreathe. Decompression required 10 minutes. Exercise stressed upper body. This 
was the first of a multiple exposure series separated by 18.0 hours at sea level conditions. 
N = 32 

16a. Same as #16 except females were tested. N = 32 

17. Same as #16 except crew returned to 14.7 psia for 18.0 hours prior to their second 6.0 hour 
exposure to 7.8 psia. No prebreathe prior to the 10 minute decompression. N = 31 
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17a. Same as #17 except females were tested. N = 31 

18. Same as #16 plus #17 except crew returned to 14.7 psia for 18.0 hours prior to their third 

6.0 hour exposure to 7.8 psia. N = 31 

18a. Same as #18 except females were tested. N = 29 

19. Exposure to 9.0 psia for 6.0 hours without prior oxygen prebreathe. Decompression required 

10 minutes. Exercise stressed upper body and 50% O 2 - 50% N 2 was used. Pilot study using 
susceptible subjects. % VGE data in these tests were reduced by 1/2 since the tested group 
was selected for its susceptibility to develop VGE. A randomly selected group would 
probably have developed the listed % VGE with the given prebreathe procedure. N = 16 

20. Exposure to 10.0 psia for 6.0 hours without prior oxygen prebreathe. Decompression re- 
quired 10 minutes. Exercise stressed upper body and 50% O 2 - 50% N 2 was used. Pilot study 
using susceptible subjects. % VGE data in these tests were reduced by 1/2 since the tested 
group was selected for its susceptibility to develop VGE. A randomly selected group would 
probably have developed the listed % VGE with the given prebreathe procedure. N = 8 

21. Exposure to 8.5 psia for 6.0 hours without prior oxygen prebreathe. Decompression required 
10 minutes. Exercise stressed upper body. Pilot study using susceptible subjects. %VGE 
data in these tests were reduced by 1/2 since the tested group was selected for its susceptibil- 
ity to develop VGE. A randomly selected group would probably have developed the listed 

% VGE with the given prebreathe procedure. N = 9 

22. Exposure to 9.5 psia for 6.0 hours using 50% O 2 - 50% N 2 mixture without prior oxygen 
prebreathe. Decompression required 10 minutes. Exercise stressed upper body. Pilot study 
using susceptible subjects. % VGE data in these tests were reduced by 1/2 since the tested 
group was selected for its susceptibility to develop VGE. A randomly selected group would 
probably have developed the listed % VGE with the given prebreathe procedure. N = 6 

22a. Same as #22 except females were tested. N = 11 

22b. Same as #22 except males were tested, (extended test) N = 20 

23. 6.0 hours oxygen prebreathe at 14.7 psia prior to 6.0 hour exposure to 4.3 psia. 
Decompression required 10 minutes. Exercise stressed upper body. N = 19 

23a. Same as #23 except females were tested. N = 19 

24. 8.0 hours oxygen prebreathe at 14.7 psia prior to 6.0 hour exposure to 4.3 psia. Decom- 
pression required 10 minutes. Exercise stressed upper body. N = 8 
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25. 2.0 hours oxygen prebreathe at 14.7 psia prior to 24.0 hours at 10.2 psia. 15 minute de- 
compression from 14.7 psia to 10.2 psia was included as a portion of the 2.0 hour oxygen 
prebreathe. 10 minute decompression from 10.2 psia to 6.0 psia after 24.0 hours. Subjects 
exercised 6.0 hours while breathing 60% 0 2 - 40% N 2 mixture. Exercise stressed upper body. 
Gas composition at 10.2 psia was 28% 0 2 - 82% N 2 . N = 15 

25a. Same as #25 except females were tested. N = 14 

26. Exposure to 8.3 psia for 6.0 hours without prior oxygen prebreathe. Decompression re- 
quired 10 minutes. Exercise stressed upper body and 50% oxygen and 50% nitrogen 
was used. N = 20 

26a. Same as #26 except females were tested. N = 11 
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